Electric dipole moment enhancement factor of thallium 
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The goal of this work is to resolve the present controversy in the value of the EDM enhancement 
factor in Tl. We have carried out several calculations by different high-precision methods, studied 
previously omitted corrections, as well as tested our methodology on other parity conserving quan- 
tities. We find the EDM enhancement factor of Tl to be equal to -573(20). This value is 20% 
larger than the recently published result of Nataraj et al. [Phys. Rev. Lett. 106, 200403 (2011)], 
but agrees very well with several earlier results. 
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A number of extensions of the standard model of par- 
ticle physics predict electric dipole moments (EDM) of 
particles that may be observable with the present state- 
of-the art experiments [l[ making EDM studies a remark- 
able tool in search for new physics. The EDMs arise 
from the violations of both parity and time-reversal in- 
variances. The present constraints on the EDMs are al- 
ready within bounds predicted by some theories [![. If 
the EDMs are not observed in the next generation of the 
experiments, some of the low-energy supersymmetry and 
other theories will be ruled out. The standard model pre- 
dicts tiny electron EDM, d e < 10 -40 e cm since it can 
not originate even from three- loop diagrams 0. 

The electron EDM is enhanced in certain atomic and 
molecular systems, and two of the most stringent limits 
on the electron EDM d e were obtained from the exper- 
iments with 205 T1: d e < 1.6 x lCT 27 e cm Q, and with 
YbF molecule: d e < 1.05 x 10~ 27 e cm Q. These limits 
significantly constrain supersymmetric and other exten- 
sions of the standard model 

Both results crucially depend on the calculated values 
of the effective electric field on the valence electron. In 
the case of Tl this effective field is proportional to the 
applied field E , E eS = KE Q , and d at ( 205 Tl) = Kd e . 
The quantity K is referred to as the EDM enhancement 
factor. 

Until recently, there was a consensus that the value of 
K is close to —580 but the latest calculation Q 

gave the value —466(10), or more than 20% smaller. All 
three of these calculations used high-accuracy methods 
that include some parts of the correlation corrections to 
all orders. Liu and Kelly [5j used relativistic coupled- 
cluster approach, but had to make various restrictions in 
their calculations to make it manageable with the com- 
puter power available in 1992. Dzuba and Flambaum Q 
used a combination of configuration interaction method 
(CI) with many-body perturbation theory (MBPT) start- 
ing from the [Xe]4/ 14 5d 10 Dirac-Fock (DF) potential and 
considering thallium as a system with three valence elec- 



trons. This potential is referred to as V N ~ 3 , where N is 
the total number of electrons. Nataraj et al. Q used rel- 
ativistic coupled-cluster (RCC) method with single, dou- 
ble and perturbative triple excitations of the Dirac-Fock 
wave functions starting from the [Xe]4/ 14 5d 10 6s 2 poten- 
tial. In this potential, referred to as V* -1 , 6s 2 shell is 
included in the core and thallium is considered to be a 
monovalent system, such as an alkali-metal atom. Both 
relativistic coupled-cluster method, in its various imple- 
mentations, and CI+MBPT method have been used for 
a number of years in many other applications, including 
study of parity violation, calculation of other (P, T-odd) 
effects, search for variation of fundamental constants, and 
many others. We note that calculations of the effective 
field for such a heavy molecule as YbF is more difficult 
than for atomic Tl, and discrepancy in the theory for 
atoms may compromise molecular limit Q as well. 

Because of the importance of this issue, we return to 
the problem of Tl EDM in this letter. We have car- 
ried out several calculations by different high-precision 
methods in different potentials. Below, we briefly sum- 
marize the calculations that we have performed and our 
main findings before providing more details of the meth- 
ods used in this work. 

1. CI+MBPT calculation in V N ~ 3 potential. 
Firstly, we have repeated the calculation of Dzuba and 
Flambaum [|| and ensured that we agree with their value 
at the same level of approximation. Then, we have cal- 
culated a number of corrections that were omitted in 
Q, including structure radiation, core-Brueckner, two- 
particle, selected three-particle, and normalization cor- 
rections. We found that some of these corrections are 
large, 5 — 7% percent, but partially canceling, causing 
the lower accuracy of V N ~ 3 results than was previously 
expected. 

2. Cl+all-order calculation in V N ~ 3 potential. 

Recently, we have developed the relativistic method Q 
combining CI with relativistic coupled-cluster (CC) ap- 
proach. This method, first suggested in (here re- 
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ferred to as the Cl+all-order method) was successfully 
applied to the calculation of divalent atom properties in 
Refs. d, This calculation allowed us to evaluate the 
effect of higher-order core-core and core-valence correc- 
tions to the EDM. We found that the effect of these cor- 
rection to the enhancement factor is small, 0.7%. 

3. CI+MBPT calculation in V w 1 potential. 
We have repeated the entire CI+MBPT calculation de- 
scribed above, including all corrections, in V N ^ poten- 
tial. We note that the CI+MBPT approach still allows us 
to fully treat all three valence electrons using CI. There- 
fore, we can accurately treat the contributions of the 
6s6p 2 states to K on the same footing as 6s 2 ns terms 
unlike the approach of Q. Our final result is based on 
the V^ -1 calculation with the higher-order corrections 
estimated from the V w_3 Cl+all-order calculation. 

4. RCC calculation of the 6s 2 ns contributions to 
EDM. In an attempt to compare with RCC calculations 
of 01) we carried out full relativistic coupled-cluster cal- 
culations including single, double, and perturbative triple 
(RCCSDpT) excitations of the DF functions in V 1 *- 1 po- 
tential. We have also used this calculation to evaluate the 
effect of the Breit interaction to EDM and found it to be 
negligible (0.36%). We have verified that our RCCSDpT 
value for the 6s 2 7s and 6s 2 8s contributions to the EDM 
are in excellent (2%) agreement with our final CI+MBPT 
contribution confirming the accuracy of our calculations. 
We have conducted RCC calculations with a truncated 
basis that we have constructed using parameters of 
The truncation of the basis to the size of Ref. pro- 
duced very large reduction (by 18%) in the size of the 
6pi/2 ~~ 7 s EDM matrix element. 

We discuss these four calculations in more details be- 
low. We start with solving DF equations Hq ifj^ = s c ip c 



where Hq is the relativistic DF Hamiltonian [8j, and 
ip c and e c are single-electron wave functions and ener- 
gies. The self-consistent calculations were performed for 
the [Is, 5d 10 ] closed core and the 6s — 8s, 6p, 7p, and 
6c? orbitals were obtained in V N ~ 3 approximation. We 
have constructed the basis set 
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consisting of 166 
orbitals (22s, 22p, 21d, 20/, 13#, and lift,). In order to 
estimate the accuracy of this basis set, we repeated some 
of the calculations with significantly larger basis set con- 
sisting of 273 orbitals (35s, Up, 28d, 27 /, 21g, and 20ft) 
and found that the differences were small and well below 
our estimated accuracy. The CI space was significantly 
increased in comparison to [f|, and included (22s, 22p, 
17d, and 16/) orbitals; higher n orbitals were allowed 
fewer number of excitations. Such CI space is effectively 
complete. All MBPT and all-order calculations were car- 
ried out with inclusion of all orbitals. 

The multiparticle relativistic equation for three va- 
lence electrons is solved within the CI framework [14 1 



to find the wave functions and the low-lying energy lev- 
els: H c ff(E n )$ n = E n Q n , with the effective Hamiltonian 
defined as H e e(E) — Hpc + ^(E)- H-pc is the Hamil- 



tonian in the frozen-core approximation and the energy- 
dependent operator T,(E) takes into account virtual core 
excitations. The £(-E) part of the effective Hamiltonian 
is constructed using the second-order perturbation the- 
ory in the CI+MBPT approach [ll[ and linearized single- 
double coupled-cluster method (LCCSD) in the Cl+all- 
order approach Q . Since the valence- valence correlations 
are very large, the CI method provides better description 
of these correlations than the perturbative approaches 
such as RCC due to possible large contributions of higher- 
order (or higher-excitation) corrections. The LCCSD 
method used here is known to describe the core-core and 
core- valence correlation very well as demonstrated by its 
great success in prediction of alkali-metal atom proper- 
ties [HI]. Therefore, combination of the CI and all-order 
LCCSD methods allows to account for all dominant cor- 
relations to all orders. 

The absolute values of the three-electron binding en- 
ergy and the energy levels of the low-lying excited states 
in respect to the ground state obtained in the pure CI, 
the CI+MBPT, and the Cl+all-order approximations are 
given in Table I of the supplementary material [l(| . We 
find that the Cl+all-order improves the accuracy of en- 
ergies and reduces the error in the ground state three- 
electron binding energy to 0.2% level. 

The atomic EDM d at of the ground state of Tl is de- 
fined as 



2E- 



|D|n)(n|g d |0) 

Ea — E n 



(1) 



where D is the electric dipole moment operator. The 
operator is given by 0: 



Hd — 2eL 



0\ Z(r) 
err 



(2) 



where d e is the EDM of the electron, Z(r) is the charge of 
the nucleus and core electrons within the sphere of radius 
r, and er are Pauli matrices. In the CI+MBPT approach, 
we construct effective valence operators for all observ- 
ables of interest 17, 13] . In this work, we need effective 



operators for the electric dipole operator D g , magnetic- 
dipole hyperfine (hfs) interaction (-ffhfs)cff, and the op- 
erator (-ff(j)eff. These operators account for the core- 
valence correlations in analogy with the effective Hamil- 
tonian. We do not perform explicit summation over 
three-particle states in our approach, but use Dalgarno- 
Lewis-Sternheimer method that involves solution of the 
inhomogeneous equation with the corresponding effective 
We include additional corrections 
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operators 

beyond random-phase approximation (RPA) in the con- 
struction of all effective operators in comparison with [f| . 
These contributions include the core-Brucckncr (<j) , two- 
particle (2P) corrections, structural radiation (SR), and 
normalization (norm) corrections. Finally, we calculated 
selected three-particle (3P) corrections to the effective 
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TABLE I: The ground state three-electron binding energy 
\E V \ (in a.u.) and the energy levels of the low-lying excited 
states in respect to the ground state (in cm -1 ) for 1/ N_1 ap- 
proximation. Results of the calculations and the differences 
from the experimental values f9, 20] (in %) are presented for 
CI and CI+MBPT approximations. 





CI 




Cl-t 


-MBPT 


Expt." 


E v 


1.9809 


4% 


2.0682 


0.2% 


2.0722 


6f>3/2 


7016 


10% 


7854 


-0.8% 


7793 


7S1/2 


24649 


7% 


26328 


0.6% 


26478 


7Pl/2 


31876 


7% 


33954 


0.6% 


34160 


7P3/2 


32834 


7% 


34974 


0.5% 


35161 


6^3/2 


33762 


7% 


36106 


0.0% 


36118 


6^5/2 


33828 


7%. 


36180 


0.1% 


36200 


8S1/2 


36549 


6% 


38693 


0.1% 


38746 



TABLE II: The magnetic-dipole hfs constants (in MHz) and 
the absolute values of the reduced matrix elements of the 
electric-dipole operator |(7j|-D||7')l ( m a - u 









Theory 


Expt. 


A (MHz) 




6P1/2 


22041 


21310.8 [21] 






7S1/2 


12395 


12297(2) [22] 






8Sl/2 


3900 


3871(1) [22] 


El (a.u.) 


(7«| 


£>||6pi/ 2 ) 


1.781 


1.81(2) [23] 




<8a| 


D||6pi/2> 


0.521 






<7*| 


D||6ps/2> 


3.393 


3.28(4) [23] 




<8s| 


0||6p 3 /2> 


0.764 





Hamiltonian [llj. We find that an accurate calculation 
of different observables in I/ N_3 potential is more com- 
plicated due to the poor convergence of the MBPT. We 
present the contributions to hfs constants A for the 8 
low-lying states in the Table II of the supplementary 
material [16j . We find that many corrections beyond 
CI+MBPT and RPA are large and partially canceling. 
As a result, an agreement between final theoretical val- 
ues and the experimental results in certain cases is not 
very good. In particular, the discrepancy between the- 
oretical and experimental values of A(7s) is at the level 
of 8%. The normalization corrections are unusually large 
(~ 6%). We calculated the normalization correction by 
approximately expressing it in terms of the derivatives 
of the MBPT corrections in respect to the energy fllj j . 
It appears that different method for treatment of this 
correction needs to be developed in the case of T/ N_3 
potential. 

We find the same problem when calculating these cor- 
rection to the EDM enhancement factor in the I/ N_3 ap- 
proximation. The CI value is -584 and the CI+MBPT, 
Cl+all-order, and RPA corrections contribute only 3, 4, 
and 3, respectively. Usually these are the most important 
corrections to the valence CI. At the CI+MBPT+RPA 
level, our result is —578 and is in a good agreement with 
the value —582(20) obtained by Dzuba in Flambaum 
using the same CI+MBPT+RPA approximation in the 



I/ N ~ 3 potential. Small difference may be due to a differ- 
ent basis set and larger CI space (including I — 3 partial 
wave) in our calculations. The corrections a, SR, 2P, 3P, 
and norm are 25, —1, —22, —2, and 36. The two-particle 
and normalization corrections are large, +4% and —6%, 
correspondingly leading to the value K = —538(46). We 
estimated the uncertainty in K based on the difference 
of the relevant hyperfine constants with experiment (8% 
for A(7s)) and the total contribution of all corrections 
beyond CI (8.6%). 

In summary, we find that the corrections beyond 
CI+all-order+RPA are large; even though they partially 
cancel each other, their total contribution is significant, 
almost 7% in I/ N ~ 3 potential. At the same time, the 
all-order CC corrections due to higher-order core- valence 
correlations are very small, 0.7%. We conclude that the 
size of different corrections to the EDM in the V N ~ 3 ap- 
proximation is not typical and missing higher-order con- 
tributions to the effective operators can be important. 
Because of that, we repeat calculations in the V N_1 ap- 
proximation. We already used this approximation in the 
calculation of the parity-nonconserving amplitude for the 
6pi/2 ~ 6^3/2 transition in Tl with 3% accuracy (24 1. 
Comparison of the I/ N_1 and V N-3 potentials for Tl cal- 
culations has been recently discussed in Ref. (25l |. 

CI+MBPT calculation in the V N_1 potential follows 
the same procedure as the one in the V N_3 approxi- 
mation, but the self-consistency DF procedure is carried 
out for the [Is, 5d 10 , 6s 2 ] core. We note that we use 
the Brillouin-Wigner variant of the MBPT in both cases. 
In this formalism, the effective Hamiltonian for the va- 
lence electrons is energy- dependent. It was shown in our 
work 26[ that the accuracy of the theory can be improved 
by calculating the Hamiltonian at the optimal valence en- 
ergy for Tl, which was found to be —1.8 a.u.. In Table |H 
we present the absolute values of the valence energy of 
the ground state and the energy levels of the low-lying 
excited states counted from the ground state obtained in 
the pure CI and in the CI+MBPT approximations. We 
note that the Cl+all-order formalism is presently lim- 
ited to the V N ~ 3 potential. In F N_1 , so-called subtrac- 
tion diagrams have to be included consistently at the 
all-order level which so far has not been implemented. 
Since the all-order core-valence corrections contributed 
only 0.7% in the T/ N_3 approximation, these are small 
at the present calculation as well. A comparison of the 
results presented in Tables U shows that the energy lev- 
els found in V N_1 approximation turn out to be closer 
to the experimental results than in V N ~ 3 approximation, 
which is already observed at the stage of pure CI ap- 
proximation. As a result, the MBPT corrections that 
give the main contribution to the uncertainty budget are 
smaller leading to better agreement between the theo- 
retical and experimental energy levels. Our values for 
the magnetic-dipole hfs constants and El transition am- 
plitudes between low-lying levels in V N_1 approximation 
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TABLE III: The breakdown of different contributions to our final value of the EDM, 
CI value and the following columns give various corrections described in the text. 


V N_1 potential. 


First column 


gives the 


CI CI+MBPT RPA Sbt 2P a SR Norm Final 


Ref. [5] 


Rcf. [6] 


Rcf. [7] 


-593.6 8.7 -13.0 16.5 -18.8 22.5 0.0 5.2 -573(20) 


-585(30-60) 


-582(20) 


-466(10) 



TABLE IV: The contributions to the EDM in our final K N_1 
CI+MBPT calculation. Columns D and Hd give reduced ma- 
trix elements of the electric-dipole and EDM operators. The 
results of our RCCSDpT calculation for the corresponding 
Qs 2 ns contributions are given in the rows labeled K(CC) for 
comparison. 



State 



AE th AE, 



6s'7s %J 2 
6s 2 8s 2 S 1/2 
6s6p 2 4 P 1/2 
6s6p 2 2 P 1/2 
6s6p 2 4 D 1/2 ? 



cxpt 



D 



H d K K(CC) 



Other 
Total 



26328 
38693 
46281 
69218 
79830 



26478 
38746 
45220 
67150 



-1.798 17.7 

-0.526 9.2 

-0.427 45.1 

-2.472 28.0 

-0.142 -4.2 



-216.6 -212.2 
-22.4 -22.8 
-74.6 

-179.2 
1.3 
-81.1 

-572.5 



are compared with experimental results 2lH23| in Ta- 
ble |TTJ The calculation of these properties was discussed 
in detail in Ref. 24j . The corrections that are taken into 
account are similar to those for the V N_3 approxima- 
tion. The only essential difference is an appearance of 
the subtraction diagrams (Sbt) in the former case. The 
differences between the theoretical and experimental re- 
sults do not exceed 3% for all relevant properties. 

In Table HTT| we present the breakdown of different con- 
tributions to the atomic EDM. The difference between 
the CI value, —594, and the final value, —573, is only 
3.7%. It demonstrates that the interaction between va- 
lence electrons is much more important and should be 
treated as accurately as possible. In Table ITVl we list the 
most significant contributions of certain valence states 
that were calculated using Eq. ([T]) and results of our 
y N_1 CI+MBPT calculation. These contributions give 
85% of the total value. Table ITVl illustrates that it is very 
important to accurately account for the contributions of 
the 6s6p 2 configurations. Their contribution to the EDM 
is ~ 45%. In the RCC method of Q, these contribution 
are accounted for as the excitations of the core electrons 
which is unlikely to provide required accuracy. 

We have also carried out completely different set of cal- 
culations using relativistic coupled-cluster method with 
single, double, and perturbative triple (RCCSDpT) exci- 
tations in y N_1 potential lU 27| to evaluate the domi- 
nant contributions to the EDM from the 6s 2 7s and 6s 2 8s 
states by a different approach. All nonlinear terms were 
included at the SD level. This method is theoretically 
very close to that of [7]. While there are differences in 
the treatment of triple excitations between the two ap- 
proaches, we find that contributions of the triple excita- 



tions to the EDM is small (less than 2%). We have also 
used this calculation to evaluate the effect of the Breit 
interaction to EDM and found it to be negligible (0.36 % 
for the 6s 2 7s contribution). Our final values for the 6s 2 7s 
and 6s 2 8s contributions are given in the last column of 
Table [TV] RCCSDpT values are in excellent (2%) agre- 
ment with our final CI+MBPT values. The agreement 
of the results obtained by two completely different ap- 
proaches confirms the accuracy of our calculations. How- 
ever, the value of this contribution inferred from Fig. 2 
of is 10% lower, about -188. We find that this dif- 
ference may be due to significant truncation of the basis 
set used in RCC calculation of 0] ■ The main part of our 
RCCSDpT calculation was carried out with very large 
numerically complete basis set (650 orbitals with I < 7). 
Fig. 2 in [3| shows that their 9s orbital already belongs to 
continuum, which is due to use of only ns orbitals up to 
14s in the RCC calculation. We have conducted a basis 
set test truncating our 166 orbital basis to n = 14 for all 
partial waves and using it in the RCC calculations. We 
find that the value of the 6pi/2 — 7s EDM matrix element 
was reduced by 18% due to basis set truncation. 

To conclude, we calculated the EDM enhancement fac- 
tor to be equal to —573(20). The uncertainty is, some- 
what conservatively, assigned based on the accuracy of 
the relevant hyperfine constants (3%) and total size of all 
corrections beyond the CI, which is 3.7%. This value dif- 
fers by 20% from the recently published result of Nataraj 
et al. 7], but agrees well with the results obtained 
by Dzuba and Flambaum Q and Liu and Kelly 
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